Large Berry phases in layered graphene 
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Brillouin zones of graphene systems possess Dirac points, where band degeneracies occur. We 
study the variety of (and large magnitude) phases that the electronic states can acquire when a 
uniform time-dependent electric field carries the electrons around one or more Dirac points in a non- 
concentric fashion. An experimentally accessible determination of excess Berry phases is proposed 
involving the Zitterbewegung of electronic current near an orthogonality point in adiabatic motion. 
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I. INTRODUCTION 

The recent research interest in graphene owes much 
to the remarkable electronic properties near the charge- 
neutrality, Dirac points (DP's, conical intersections or 
linear electronic degeneracies) at the K and K'-points in 
the Brillouin zonei^^^. Also several studies devoted 
to graphene layers of various thicknesses have appeared 
recently Z&tL i 10 : 11 : 12 : 13 . Based on the earlier elucidation 
of the Brillouin-zone structure in Refs. [ll[ll and the 
more recent works in Refs. Il6lfl9l in bilayered graphene 
one centrally and three trigonally arranged DP's around 
the K-type points can be located. The electronic mo- 
tion is described by a Hamiltonian that includes interac- 
tions between different sites in either plane and between 
the two planes and expressed the coupling between four 
electronic bands. An approximate and simplified model 
using two 2x2 Hamiltonians (for the K and K'-points, 
respectively) was given in Ref. [H. The repercussions of 
different Hamiltonians on the zero energy minimal con- 
ductivity were discussed in Ref. [HI, where the effect of the 
trigonal distortion was found to be significant, whereas 
the simplified Hamiltonian of Ref. [l^ led to results similar 
to those of the four band model. 

Another development regarded the Berry phas e 20 : 21 
acquired upon circling around the DP's in the k- 
plane&i&H. It was assumed in these works that a cir- 
cling around each point leads to an added value of ±7r to 
the geometric phase, with the sign depending on deriva- 
tives in the neighborhood of the intersection points. (The 
sign of the Berry phase encountered in the molecular 
physics, electron-vibration context was determined in 
Refs. [22l l 23l ) It is the purpose of the present work to show 
that when a Berry phase is actually created in a physical, 
adiabatically cyclic process, then the magnitude of the 
phase will differ from the assumed value of ±7r and is 
more correctly (2n + 1)tt, where n is a signed integer or 
zero. The above result holds for a single DP, when the 
circling is performed in a non-concentric manner (and, 
likewise, when the circling is around an odd number of 
intersections). Analogously, when the circle is around 



two or an even number of intersection points the assumed 
result of ±27r or zero is superseded by 2n7r (with n tak- 
ing values, as before). These predictions are based on 
a proper treatment of the time- dependent adiabatic pro- 
cess, which requires near an orthogonality point a formal- 
ism that involves a correction term beyond the extreme 
adiabatic limi t 24 ' 25 : 26 . (At an "orthogonality point" the 
momentary wave function has zero overlap with the start- 
ing wave function. The integer n is related to the vary- 
ing speed of the electron along its cyclic path. This is 
quantified by the parameter n', named "the inverse rel- 
ative slowness" and defined below in equation (Til?)) . In 
the present context and formalism the extreme adiabatic 
limit is defined by the vanishing of the ratio e in equa- 
tion (14) for all times, while n' in equation (18) is finite 
and possibly numerically large. A conventional treatment 
would first neglect the second and third terms in equa- 
tion (16) and then obtain a ±7r Berry phase from the sign 
change of the first term. We retain the remaining terms 
and calculate the acquired phases at the several zeros 
near the orthogonality point, neglecting the small terms 
only after exiting from the orthogonality neighborhood.) 

Now it is true that Berry phases are usually written 
modulo 2-7T, and are gauge invariant only under conditions 
of complete cycling, but the effect of the extra \n\ 7^ 
shows up even before completion of the cycling, as an 
open path phased!, and can in principle be measurable. 
(The extra acquired phase is gauge invarian t 24 : 25 : 26 . Ex- 
perimental verification of the extra open path phase has 
been proposed before^.) An essential point (not present 
in previous electron-nuclear settings) is that in graphene 
it is technically feasible to control the phase acquisition 
by varying the electron concentration and by an elec- 
tric field. The former is achieved by the manipulation of 
the (time-constant) gate voltage (through the substrate) 
and the latter by a well designed application of a time- 
varying uniform in-plane electric field that guides the 
electronic motion through the " acceleration theorem" (of 
which more in the sequel) . Application of a time indepen- 
dent electric field was treated in Ref. Hi| In the present 
setting of the adiabatic change in the electronic wave vec- 
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tor, the magnitude of the applied electric field would be 
of the order of 10 4 V/m or less; this is rather less than 
the value 10 7 V/m proposed in the preceding reference, 
appropriate to conditions that are not adiabatic. 

Although in this work we do not detail the experimen- 
tal aspects of the phase acquisition, we do propose and 
investigate theoretically an entirely new method for the 
experimental observation of the Berry phase, feasible for 
high Berry phase values. (In our future mentioning of the 
phase, the dynamic phase of Ref. is understood to be 
subtracted from the total wave-function phase. We re- 
turn to the dynamic phase in section fVIin This is made 
possible during the acquisition of the phase in the vicinity 
of the orthogonality point and is measurable through the 
Zitterbewegung (ZB) of the electronic motion (current). 
ZB in graphene related materials was the subject of sev- 
eral work o 33 i 34 ' 35 i 36 . (Also related is referenced?].) It is 
shown in section IVll that in a adiabatically time develop- 
ing eigenstate, ZB is observable near the orthogonality 
points. The time duration of the ZB is associated with 
the inverse relative slowness parameter. In the present 
context where adiabatic changes are considered, the du- 
ration for which ZB has to be observed is expected to 
be rather large, of the order of 10 3 fs or more. (Restric- 
tions on the duration of the ZB by the finite width of 
the electronic wave-packet were previously considered in 
Rcfs. 34.36.38.) 



II. GRAPHENE-BASED HAMILTONIANS 

A. Monolayer graphene 

This is a zero-gap semiconductor, whose zero- 
temperature electronic properties in the undoped or 
slightly doped form arise from the neighborhood of 
hexagonally arranged K-points (the Dirac points, DP) 5 . 
Here the valence and conduction bands meet in a coni- 
cal intersection. The Hamiltonian, expressed by means 
of the planar quasi-particlc momenta (k x ,k y ) measured 
from a K-point, have the form of a two-dimensional 
Dirac-Hamiltonian 



H = v F 



- k i- 



(1) 



where a real representation was adopted for the electronic 
band states and vf is the Fermi velocity (about 10 6 m/s). 
The off-diagonal, complex Hamiltonian commonly used 
in the graphene literature is M'HM, having applied the 
unitary transformation matrix 



M 



(2) 



to the above Hamiltonian H. 

For future use we note that for momenta \k\ whose 
cyclic time rate of change is represented by to, the require- 
ment of adiabaticity is given in a monolayer graphene by 

V F \k\^>LO. (3) 



B. Graphene bilayer 

Four electronic bands describe the salient properties. 
A concise form of the 4-band Hamiltonian for an electron 
in a k (wave vector)state near the K and K'-points was 
written out in Ref . [ijj (Eq. (1) there.) An approximate 
2-band Hamiltonian proposed by Ref. [H can be written 
in a real electronic state representation and following the 
real-matrix notation used in Refs. I24lf26l as 



-U(k) V(k) \ 
V(k) U(k)J 



(4) 



for (e.g.) the K-point. Here the designation k is used 
for the reduced in-plane wave-vector measured from the 
K-point. 



M) 



k (k x , k y ) _ ~ - 



k 



(5) 



where ko 



27173 



having introduced the intra-layer cou- 



pling strength 70 between the basis atoms A and B within 
the first layer and between A' and B' within the second 
layer, the vertical interlayer A-B' coupling strength 71, 
the (weaker) diagonal interlayer A'-B coupling strength 
73 and the honeycomb lattice constant a. Numerical val- 
ues that have been accepted for these parameters are 
7o=3.16 eV, 7i=0.39 eV, 7 3 =0.315 eV, 0= 0.246nm, fc = 
0.05775 ran -1 . The matrix elements are 



U(k) = eo [k x -(k 2 x -k% (6) 
V(k) = e [k y + 2k x k y ], (7) 

with an overall coupling energy cq = 71 (^) 2 = 3.87 

mcV. For a K'-point the signs of the first term in U (k) 
and of the second term in V(k) need to be changed. 



III. MOTION OF ELECTRONS 

When subject to an electric field E, the wave vector of 
an electron (with charge — e) changes in time (t) accord- 
ing to the "acceleration theorem" 30 ' 31 



dk 



-E. 



(8) 



A spatially uniform and time varying electric field is ex- 
pressed in terms of a vector potential A(t) through 



E 



dA(t) 
dt 



giving 



k(t) = -A(i) 



(9) 



(10) 



apart from an initial value of the wave vector. Thus, 
the above Hamiltonians are to be understood as func- 
tions of the vector potential or of the time integral of the 
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externally applied electric field. This will be the mean- 
ing attached to the independent variables (k or k) in 
the Hamiltonians, though (for simplicity of notation) we 
shall continue to write them as functions of the time- 
dependent planar wave-vectors, rather than that of A. 
In the sequel we shall assume the following time depen- 
dence for the moving reduced wave vector: 

k x (t) = acos(uit)+c, (11) 
k y (t) = 6sin(w£), (12) 

with constants a, b and c and a period of T = 2tt/uj. 



IV. GEOMETRIC PHASE OF THE WAVE 
FUNCTION COMPONENT 



Solutions of the time dependent Schrodinger equation, 
written for a general form of a periodic 2x2 Hamilto- 
nian as in equation Q and valid in the adiabatic limit, 
were given in Refs. [24l26l . with special regard to the ge- 
ometric phase acquired near an " orthogonality time" t v . 
This is the instant at which the initially engaged compo- 
nent of the quasi-spinor vanishes, or the state becomes 
orthogonal to the initial state. There may be several 
such instants, depending on the path taken by the time 
dependent parameters. As already remarked in the In- 
troduction, the solutions of the above cited works went 
one approximation (in the adiabatic parameter, to be 
defined shortly) beyond the formal adiabatic wave func- 
tion. The instantaneous energy and twice the mixing an- 
gle appropriate to the Hamiltonian in equation ([4| (and 
analogously to that in equation |T]) ) are 



W(t) = ^£/(k(i)) 2 + W*)) 2 , (13) 
V(k(t)) 



x(t) = arctan ■ 



(14) 



U(k(t)) 

The (small) adiabaticity parameter is in terms of these 

*(*) 



e(t) 



w(ty 



(15) 



We now turn to the solution of the time dependent 
Schrodinger equation in the form 



£ ( /(*) 

dt{g(t) 



-U(t) V(t) 
V(t) U(t) 



(fit) 

\ g(t) 



(16) 



with the initial condition for an energy eigenstate f(t = 
0) = 1, g(t = 0) = 0. In Refs. HHl we have obtained 
(after the removal of the dynamic phase factor) the fol- 
lowing component amplitudes of the wave function in the 



neighborhood of an orthogonality point t v as 



1 



2\W(t v )\ 



(l-Tm'^e-W^ie-'"))] 



1 



g{t) = - sm-x{t v ) = -1. 



(17) 
(18) 



The expressions are correct to the order of \x(t v )(t—t v )\ 2 

-i e(0) 



n'(t v ) = 7r 



(19) 



which ratio compares the slownesses (or adiabaticitics, 
defined in equation (|15j) ) at the starting points and at 
the orthogonality point. It is a key quantity in the phase 
acquisition phenomenon, as we shall describe below. In 
contrast, the last term on the right of equation (fT7| is 
of little importance. (This has been confirmed both the- 
oretically and by numerical computation). In it <f>(t v ) is 
the fractional part of the energy integral 2 J Q V W(t)dt. 
(See Eq. (31) in referenced.) The energy (or frequency) 
2| W(£n) | in the complex exponential is the energy separa- 
tion between the two eigenstates near the orthogonality 
point and will show up in our treatment of the electronic 
motion (or current) in section IVT1 as a Zitterbewegung. 

We shall now describe the geometric phase acquisition 
process in terms of the minority component f(t) in equa- 
tion l|17p. Were it only the first term in f(t), a phase of 
— 7r would accrue upon traversing the zero at t = t v from 
t < t v to t > t v . This phase is continuously carried in f(t) 
and would remain with the wave function also when a full 
revolution is made at which instant the g(t) component 
becomes again zero. Thus, — tt is indeed the final result 
(the Berry phase) as long as n'(t v ) is smaller than or of 
the order of 1. However, when n'(t v ) is numerically large, 
there are further zeros in the complex t-plane, which are 
all located in the lower half of the complex t-plane (at 
shallow depths of the order l/\W(t v )\)ZL&22.. Looking 
now for the moment at the behavior of the real part of 
f(t) [Re/(i)], with algebraically increasing real t, as this 
passes each adjacent (complex) zero, Re/(i) acquires a 
further — it phase. This is also how the full function f(t) 
acquires its phase, with the role of the imaginary part be- 
ing to ensure that the phase of f(t) changes in a smooth, 
continuous way. The rate of tt- acquisitions is clearly 
2|lT / (t„)|/7r (a constant for a given orthogonality point 

tv). 

How many (complex) zeros are there in the asymptotic 
limit of large \n'(t v )\? Taking into account that the mod- 
ulus of the complex exponential in equation p7[) is unity, 
we clearly see that the complex zeros arise only as long as 
\t— t v \ < J[vy(t"j| 1 smce otherwise the first term will dom- 
inate the function and zeros are not possible. Taking into 
account the previously obtained value |2W(i„)|/7r for the 
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rate of acquisition of 7r's, we see that the number n(t v ) 
(an integer!) of acquired n phases is for either positive or 
negative t—ty's about 2n'(t v ). [The same result would be 
obtained by counting the number of loops around the ori- 
gin of f(t), since each (anti-clockwise) loop contributes a 
2-7T phase and each loop around the origin implicates two 
zeros of Re f(t) (as well as of Im /(£)).] 

The net result is then that the acquired Berry phase 
around t = t v is 

[1 + 2n(t v )]ir (20) 

where n(t v ) is a signed integer (or zero) close to the ra- 
tio n'(t v ) defined in equation (|15[) . n' is named the in- 
verse slowness ratio and examples taken from molecular 
degeneracies have shown that it is unity for concentric 
circular motion around a single intersection point (and, 
in the case of several intersection points, for concentric 
circular motion around one intersection point provided 
the motion is sufficiently far away from other intersec- 
tions). However it can be large for circling that starts 
close to an intersection point and is not concentric with 
i t 24 i 25 i 26 i 32 . For graphene layers similar results for the 
phase have now been found and are shown in a set of fig- 
ures. These are obtained by solving the time dependent 
Schrodingcr equation and in them the acquired phases 
agree very closely with (2nf + 1)-7T, in which nf is given 
by the formula in equation (|19p . This depends only on 
quantities contained in the Hamiltonian. The case of 
(monolayer) graphene, with essentially isolated DP's (the 
K and K' points), is formally identical to the single de- 
generacy results shown in section 3 of reference |32| and 
will not be reproduced here although monolayer graphene 
may be the more convenient candidate for the experimen- 
tal verification of the theory (However, the one order of 
magnitude lower value of W(t v ) in bilayer graphene than 
in monolayer may make the oscillations easier to moni- 
tor.) 

In a graphene bilayer we have a DP at each K (and 
K') with three satellite DP's, trigonally situated around 
each if- type point. Thus, this is a new situation; also 
for the reason that there exist for it both two and four 
dimensional Hilbert space description s] 10 ' 12 . The Hamil- 
tonian in the former description is shown in equation (J6]) 
and equation (J7J) ; some essential features of the energy 
contours are shown in Fig. [T] 

We begin by showing results for a contour that starts 
close to and just outside a trigonal DP and makes a bee- 
line with a path ultimately skirting round the central 
DP. [In Koshino and Ando's reference [1(3 figure 2, the 
path might start at an energy of about (or smaller than) 
0.5eo (roughly equivalent to an electron density below 
4xl0 10 cm~ 2 ) or, otherwise, at k x < 1.2fco, k y = and 
then return so as to come down beyond the origin of 
k, the K-point.] A full contour of this type is not ex- 
pected to change the sign of the wave-function (with the 
dynamic phase disregarded), since two DP's are encom- 
passed by it. In this case one obtains the variation of 
the 2n part of the phase near an orthogonality point, 



k y 







1 

'••[ 


. 5 

1 














1 'o 


.5 






















-1 . 5 


-1 


-0|.5 
j 




0.5 


I 


1.5 






— k 


. 5 
















-i 














1 -1 


.5 











FIG. 1: (color online). Contours for bilayer graphene (in the 
two state-representation). Broken curves (in red) : The po- 
tential V — (in the Hamiltonian of equation Q). Dotted 
curve (in blue) : The potential U = 0. The circle (in green) 
denotes an evolution path. The positions of the Dirac points 
are at the crossings of the broken and dotted curves. In the 
text, the evolution starts at a positive value of the abscissa. 
Orthogonality points t v , where x = Tj are where the path in- 
tersects the broken (red) curve. For the evolution path shown 
(circling radius =0.5, shift =0.3), there is only a single orthog- 
onal position t v on the circular path, which is located on the 
negative abscissa. For circles shifted to the left it is possible 
to have two more t v 's. 



which is of interest. In figure [5] we show (with open cir- 
cle symbols) the resulting n? in a rather simple circling 
situation, when the contour starts to the right of the cen- 
tral DP and cycles anti-clockwise round it and inside the 
trigonal DP's. An orthogonality point t v is met at about 
(but not at exactly) half a full circle. Here the (open 
path) geometric phase makes a jump, whose magnitude 
is close to the ratio n'(t v ) defined above. For a circular 
contour around the central DP (which is the path here 
taken) the magnitude of nf depends on the distance of 
the starting point (located between the rightmost and 
the central DP's). As figure [2] shows, \n'\ is large for 
contours starting close to the central DP and decreases 
as the starting point recedes only to become again large 
when the starting point is precisely at half way to a trig- 
onal DP. This finding arises from the presence of several 
DP and is absent for a Brillouin zone with a single DP. 
The triangles in the same figure show, for comparison, 
the phase jumps as obtained by numerically solving the 
time dependent Schrodingcr equation (the subject of the 
next subsection). 
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FIG. 2: (color online). Inverse relative slowness parameter n' 
for adiabatic circling plotted by open circles (in blue) against 
the starting distance from the central Dirac point (DP) in 
the direction of a trigonal DP (given in units [DP] , measuring 
the distance of the trigonal DP from the central one), n is 
approximated by the integer n in the Berry phase (2n+l)7r ac- 
quired at an orthogonality point t v during circling. A circular 
path with radius 0.5[DP] was assumed. The inverse slowness 
parameter n' is large when the path starts near the central 
DP or (numerically) at just half way to the trigonal DP. (Here 
n' turns negative.) The triangles (in red) show the results of 
numerically computing the phase jumps. Some triangles are 
missing because of numerical difficulties in achieving conver- 
gence. 



NUMERICAL VERIFICATION OF THE 
"INVERSE RELATIVE SLOWNESS" 
PREDICTION 



V. 



In Refs. l24H26l the Berry-phase was associated (through 
what we regard as the proper treatment of the time de- 
pendent Schrodinger equation (TDSE) in the near adi- 
abatic limit) with the quantity (2n' + 1)7T (or 2n'ir), n' 
being (as described above) the inverse relative slowness. 
In the quoted articles the theory was verified by numeri- 
cal solution of the TDSE in the adiabatic limit, in the 
context of nuclear motion in electronically degenerate 
molecules. In the present graphene context (which differs 
from the molecular setting in essential physical, and some 
formal, details) the numerically obtained phase jump val- 
ues n have been plotted in Fig. [2] with triangles. These 
show good agreement with the values (open circles) ob- 
tained algebraically. Where the triangles are missing, 
this is due to numerical difficulties in achieving conver- 
gence. 

An analogous verification has been performed by solv- 
ing the TDSE for a graphene bilayer Hamiltonian, in- 
cluding a periodic, time dependent electric field E(t) and 
calculating at each instant t the phase of the wave func- 
tion. In figure [3J the phase development is shown for 
one full cycle of the electric field and the detail for the 




FIG. 3: Adiabatic phase evolution in bilayered graphene 
(obtained from numerical solution of the time dependent 
Schrodinger equation within the two-state formalism). The 
adiabatic revolution period T is 27r x 10 5 in inverse units of 
the Fermi velocity times the DP distance. Circular cycling 
meets three orthogonality points along its path (as explained 
in the caption of Fig. [1]). In the upper drawing, at these points 
steep rises take place in the phase, with the middle one (on 
the left-part of the abscissa in figure [1} being the longest rise. 
The near-horizontal wavy part becomes straightened out, as 
the motion becomes more adiabatic. In the lower drawing the 
first t v region (above the k x axis) is enlarged. 



first step. The non-constant, rising parts of the phase 
near the orthogonality points will be the subject of the 
section on " Currents" . 

The drawings in figures [T][3] have been based on the 
two-state (approximate) description of bilayer graphene. 
In figure [4] we compare the component phases ob- 
tained from forward integration of the time dependent 
Schrodinger equation for the four-state with that for the 
approximated two-state description^. There are no dif- 
ferences visible to the naked eye, which we interpret as a 
satisfactory test for the robustness of the adiabatic the- 
ory formulated in our previous papers. 



6 




I ' 1 1 1 1 1 1 1 1 

0.0 0.25 0.5 0.75 1.0 

t/T 









ra 


150 










o 












a 


100 










a 






o 






• i-H 

Sh 


50 
















£ 







o 




0> 






fed) 







-50 1 1 1 1 1 1 ' 1 

0.0 0.25 0.5 0.75 1.0 

t/T 

FIG. 4: (color online). Comparison of the Berry phase be- 
haviors of the four- and two-state models in graphene bilayer. 
Curves for the four state model are slightly displaced down- 
ward, for clarity. Two different starting positions to the right 
of the central Dirac point (DP) were chosen, expressed in units 
of the distance to the right-most DP . Upper frame starting 
position: 0.8. Lower frame starting position: 0.25. The adi- 
abatic revolution period T is 2tt x 10 4 in inverse units of the 
Fermi velocity times the DP distance. Note the different or- 
dinate scales. 



VI. ZITTERBEWEGUNG (ZB) OF THE 
CURRENT 

The equation of motion for the Heisenberg position op- 
erator r(i) of an electron is ir(t) = —[H,r(t)]. The cur- 
rent is j(t) = — er(t). We shall evaluate this for a mono- 
layer graphene near its K point where the expressions 
are quite simple, rather than for a graphene bilayer, for 
which the results are analogous, but more complicated. 
Inserting the Hamiltonian from equation ((T|) one obtains 
immediately the velocity components as 

r(t) = x(t)l + y(t)i = v F Q (21) 



The observable, expectation value of the velocity is the 
expectation value over the wave function, namely, 

<*(t) <f)Q i) (£)' ( 22 ) 

where the stars denote complex conjugates. Two situa- 
tions need now to be considered. 



A. Electronic motion during most of the contour 
path 

1. Extreme adiabatic limit, e —> 

Here 

(A .9*) - (/,<?) = (cosi X (*),-sinix(i)) (23) 

in terms of the mixing angle hx(t), which effects the di- 
agonalization of the 2x2 Hamiltonian. For the mono- 
layer Hamiltonian in equation |T]) this is defined through 
x(t) = arctan-^ = — arctan ^ x |*j , so that it is possi- 
ble to express the momentum components as k x (t) = 
|k|cosx(£), k y (t) = —\k\smx(t), where the time depen- 
dence arises from the applied vector potential as in equa- 
tion (fT0|) . Then equation (|22|) evaluates to 

< r(i) >= v F [cos X (t)i- sinx(t)j], (24) 

which merely reaffirms the equivalence of the time be- 
haviors of the momentum and velocity vectors for a com- 
pletely adiabatic motion along most of the pathway. 

2. First order adiabatic correction 

When the wave function components are expanded cor- 
rect to the order of the small adiabatic parameters e(i), 
one obtains the terms describing the ZB of the electronic 
velocities. 

<i(t)>ZB =-^lsm X (t)sm2I(t), 
<y{t)>ZB =-^lcosx(i)sin2J(i), (25) 
where we have denoted the integral over energy as 

I(t) = f dt'W(t'). (26) 
Jo 

There may be some practical difficulties in the obser- 
vation of this motion as it is superimposed on a larger 
motion shown in equation l]24l) . 

B. Near the orthogonality point t v 

Here the eigenstate component amplitudes (/, g) take 
the time-dependent forms shown in equation (jTTJ) and 
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equation (jl8|) . Substitution into equation (|22|) shows 
that in the adiabatic limit the x-component moves with 
a uniform speed of Vf- However, the y-component of the 
electron has the value 



y{t) = ±2v f Ref(t). 



(27) 



As discussed in section IIV1 this is a small quantity, of 
the order of the small adiabaticity parameter e(t v ) which 
oscillates with the period characteristic of ZB, given in 
the present case by ir/W(t v ). The eigenstate oscillations 
take place for a time interval of about ^ ^ , placed sym- 
metrically about t v . This interval can be quite large in 
case of strongly non-concentric circling, for which the rel- 
ative inverse slowness n' (defined in equation (|19p ) has 
large values (e.g., of the order of 10 2 )^. 



VII. SUBTRACTION OF THE DYNAMIC 
PHASE 

Results in the previous sections have the following 
implication (again phrased for simplicity's sake for the 
monolayer graphene case, section Til A[) : For a given ra- 
dius of circling k c around a K-point the acquired Berry 
phase can be increased by diminishing the distance Ak 
between the starting point and the center of the circular 
path. The increase in the Berry phase is about 2tt times 
n' (the relative inverse slowness) which is unity when 
the distance Ak — k c and increases beyond bounds as 
Afc — > 0. A hypothetical direct experimental determina- 
tion of the augmented phase would then likely perform 
two measurements at or close to these two limits. 

Since the above theory relates to the Berry phase, 
while the total acquired phase includes also the dynam- 
ical phase (designated (/>d), it is purpos eful to estimate 
the latter. It was shown in Rcfs. I24H26I that in the adi- 
abatic limit [equation ([3]) above] this is simply given by 
the time integral of the instantaneous energy. Then for 
a uniform circular motion in the k-plane 

4> D {k c ,Ak) 

/>2tt/lo 

= v F / V( k c cosojt) - Ak) 2 + (k c sin u;t) 2 dt. (28) 
Jo 

This elementary integral has the following values in the 
two limiting cases discussed above: 



concentric cycling: <pD{k c , k c ) = 2tt 



v F \k c 



touching: </>£)(fc c ,0) = 2\/2 



LU 

v F \k c 



(29) 
(30) 



These dynamic phase values can be subtracted from the 
observed total phase to obtain the Berry-phase. From 
our direct computations of the acquired total phase by 
the wave function, for situations such as shown in Fig. [3] 
we find that for large \n'\ -values the phase is dominated 
by the Berry-phase. (The dynamic phases at the three 
phase jumps shown in that figure amount to -20, -60, -20 
radians, respectively.) 

VIII. CONCLUSION 



This work has focused on two (interrelated) issues in 
the context of mono- and bilayered graphene. 

First, (based on our past works on adiabatic cycling, 
which are only briefly recapitulated here) we have shown 
that large amplitude (and gauge-invariant) phases of 
the electronic wave- functions can be generated in two- 
dimensional structures of the graphene-type. A concep- 
tual experimental procedure involving a uniform time 
varying electric field has been outlined, though without 
our proposing detailed prescription. The key geometrical 
element is the non-concentric adiabatic path in k-spacc 
around the K points. 

The second issue is the association of the developmen- 
tal stage of the large amplitude Berry phase (near the or- 
thogonality points) with the Zitterbewegung of the elec- 
tronic motion and its observational possibility. So far 
little attention has been given to the possibility of exper- 
imentally detecting the Berry phase in the neighborhood 
of an orthogonality point; though some theoretically ori- 
ented remarks regarding the abrupt nature of this devel- 
opmental stage of the Berry phase have been mad o 40 ' 41 . 
These remarks addressed the minimal ±7r phase acquisi- 
tions near each orthogonality point, whereas our proposal 
for the observation of the phase by Zitterbewegung is for 
large Berry phases. 
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